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A metal-insulator transition (MIT) in BiFeOa under pressure was investigated by a method 
combining Generalized Gradient Corrected Local Density Approximation with Dynamical Mean- 
Field Theory (GGA+DMFT). Our paramagnetic calculations are found to be in agreement with 
experimental phase diagram: Magnetic and spectral properties of BiFe03 at ambient and high 
pressures were calculated for three experimental crystal structures R3c, Pbnm and Pm3m. At 
ambient pressure in the R3c phase, an insulating gap of 1.2 eV was obtained in good agreement 
with its experimental value. Both R3c and Pbnm phases have a metal-insulator transition that 
occurs simultaneously with a high-spin (HS) to low-spin (LS) transition. The critical pressure 
for the Pbnm phase is 25-33 GPa that agrees well with the experimental observations. The high 
pressure and temperature Pm3m phase exhibits a metallic behavior observed experimentally as 
well as in our calculations in the whole range of considered pressures and undergoes to the LS state 
at 33 GPa where a Pbnm to Pm3m transition is experimentally observed. The antiferromagnetic 
GGA+DMFT calculations carried out for the Pbnm structure result in simultaneous MIT and 
HS-LS transitions at a critical pressure of 43 GPa in agreement with the experimental data. 

PACS numbers: 71.30.+h, 71.27.+a, 71.20.-b 


I. INTRODUCTION 

Multiferroics are used in various applications for en¬ 
ergy production, transmission of high voltage lines, data 
storage devices, and sensors— They will help to replace a 
number of currently used lead-based materials which con¬ 
tain lead toxic and harmful to the environment. One of 
the most promising candidates for applications, bismuth 
ferrite BiFeOa, is actively studied because of the cou¬ 
pling between ferroelectric and magnetic order around 
room temperature in this compound. Recent investiga¬ 
tions of BiFeOa at high pressures up to 60 GPa reveal 
metal-insulator— and high-spin to low-spin (HS-LS) tran¬ 
sitions^ in Fe 3+ at room temperature in a relatively wide 
pressure range 40-55 GPa accomplished by the structural 
phase transition— 

The problem becomes even more complicated taking 
into account that BiFeOa has a rich phase diagram. At 
ambient pressure and up to 1100 K it has rhombohe- 
dral (l?3c) crystal structural. Increasing temperature 
and pressure, the structure of BiFeOa transforms into 
orthorhombic Pbnm*2, and then cubic Pm3mi— 

For many years, metal-insulator transition (MIT) in d 
or / metal compounds^ is one of the central issues in con¬ 
densed matter physics. The most spectacular examples 
are pressure-driven transitions from a wide gap Mott in¬ 
sulator to metallic state in transition metal oxides. For 
MnO and Fe20a (d 5 configuration), the metal-insulator 
transition is accompanied by the high-spin to low-spin 
(HS-LS) transition. 

Recently, these and other MITs and spin transitions 
were successfully described theoretically^— employing 


the method combining density functional approximations 
(like, GGA or LDA) with dynamical mean-field theory— 
Pressure-driven MIT correlated with magnetic collapse 
could be treated as a delocalization of magnetic electrons 
or structural phase transition into new phase with Neel 
point below room temperature— 

Another Mott-type mechanism controlled by dramatic 
U eff decrease due to HS-LS transition was proposed in 
Ref—. Mott-type MIT driven by a broadening of the t 2g 
states was confirmed by the LDA+C7 calculations which 
revealed the HS-LS transition at 36 GPa and a transition 
to the metallic phase with no localized moment above 72 
GPa—. 

In this paper we investigate the properties of BiFeOa 
under pressure employing the method accounting for dy¬ 
namical electronic correlations and considering three ex¬ 
perimental crystal structures. 


II. METHOD 

The GGA+DMFT method^ is realized in a computa¬ 
tional scheme constructed in the following way: first, a 
Hamiltonian Hqga is produced using converged GGA re¬ 
sults for a compound under investigation, then the many- 
body Hamiltonian is set up, and finally the correspond¬ 
ing self-consistent DMFT equations are solved. In this 
work the Hamiltonians Hgga are constructed in a Wan- 
nier function (WF) basi s 18 ! 19 using the projection proce¬ 
dure described in detail in Ref. 12(1 . Initial ab-initio cal¬ 
culations of the electronic structure are done within the 
pseudo-potential plane-wave method, as implemented in 
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Quantum ESPRESSO^. 

The WFs are defined by the choice of Bloch functions 
Hilbert space and by a set of trial localized orbitals that 
will be projected on these Bloch functions. The basis set 
includes all bands that are formed by O 2 p and Fe 3d 
states and correspondingly the full set of the O 2 p and 
Fe 3d atomic orbitals to be projected on Bloch functions 
for these bands. That would correspond to the extended 
model where in addition to the d orbitals all p orbitals 
are also included. 

The resulting p—d Hamiltonian to be solved by DMFT 
has the form 

H = Hgga - H dc + i £ utftiLffa, (1) 

i,a,/3,cr,a' 

where is the Coulomb interaction matrix, nf aa is 

the occupation number operator for the d electrons with 
orbitals a or /? and spin indices a or a’ on the i-th site. 
The term Hdc stands for the d-d interaction already ac¬ 
counted for in GGA, so called double-counting correc¬ 
tion. In the present calculation the double-counting was 
chosen in the following form Hd c = U(ndm(t — \)I- Here 
fidmft is the self-consistent total number of d electrons ob¬ 
tained within GGA+DMFT, U is the average Coulomb 
parameter for the d states, and / is unit operator. 

The elements of matrix are parameterized by U 
and Jh according to the procedure described in21. The 
values of Coulomb repulsion parameter U and Hund ex¬ 
change parameter Jh were calculated by the constrained 
LDA method^ on Wannier functions^. The values U =6 
eV and Jh= 0.93 eV obtained in these calculations are 
close to the previous estimation s 16 ! 23 . The effective im¬ 
purity problem in DMFT was solved by the hybridiza¬ 
tion expansion Continuous-Time Quantum Monte-Carlo 
method (CT-QMC)2i. The calculations for all volumes 
were performed in the paramagnetic state for inverse 
temperature /3 = 1/T=15 eV -1 corresponding to 770 K. 
According to the phase diagram, no long-range ordering 
in BiFeOa was found in this temperature region. The 
spectral functions on real energies were calculated em¬ 
ploying Maximum Entropy Method (MEM)^. 

III. RESULTS AND DISCUSSION 

In all structures under investigation iron is surrounded 
by oxygen ions forming an octahedron. In the high pres¬ 
sure cubic Pm3m phase, the Fe d band is split by crystal 
field into threefold degenerate t^g and twofold degener¬ 
ate e g sub-bands. A trigonal distortion in the R3c struc¬ 
ture lowers the point symmetry group Oh to D^h, so that 
threefold degenerate t 2 g becomes split into twofold degen¬ 
erate e g and non-degenerate a\ g band. The Pbnm struc¬ 
ture has the lowest symmetry and all orbitals become 
nonequivalent but still three groups of orbitals could be 
considered: e gl and e g2 , e g i and e g2 and ai g . For sim¬ 
plicity we will use cubic t 2 g and e g orbitals notations for 
analysis hereafter. 
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FIG. 1: (Color online) Spectral function of Fe d states for 
different volumes related to the ambient pressure volume Vo 
obtained in the GGA+DMFT (CT-QMC) calculations for the 
low pressure R3c phase at 770 K. 


Previous calculations demonstrated that LDA fails to 
describe an insulating character of the BiFeOa ground 
state at ambient pressure and for all volumes it is metal¬ 
lic in all structures. However, a gap of the AFM origin 
was obtained in the GGA calculation made for several 
magnetic structures^. 

Our GGA+DMFT calculations produced for the R3c 
phase of BiFeOa show that taking into account Coulomb 
correlation effects results a wide-gap Mott insulator and 
high-spin state for ambient pressure in agreement with 
the experimental data. The obtained spectral functions 
for different cell volumes are shown in Fig. [Q The cal¬ 
culated energy gap value of about 1.2 eV agrees well 
with the experimental value of 2.4 eV for direct optical 
ga p £ £ i'— and 1.8 eV indirect optical gap at AP— at room 
temperature and indirect gap value 1.3 eV measured in 
thin films^i, 1.5 eV at 820K—. The magnitude of mag¬ 
netic moment yj< rrii > is 4.6/is at AP. This number 
agrees very well with the high-spin state of the Fe +2 ion 
(d 5 configuration) in cubic crystal field: 2 electrons in the 
e g band and 3 electrons in the t 2 g band with the magnetic 
moment value of 5 ps- Then the magnetic moment de¬ 
creases and at 0.76Vo drops down to 2.5 pBi and it equals 
to 2.0 pb when the volume is 0.72Vo (63 GPa). The later 
value is close to the one expected for the LS state with 
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FIG. 2: (Color online) Spectral function of Fe-d states vs. 
volume reduction obtained in GGA+DMFT (CT-QMC) cal¬ 
culations for the Pbnm phase at 770K. 


5 electrons in the t 2 g band but strong hybridization be¬ 
tween iron and oxygen leads to the non zero occupation 
of the e g orbitals. This mechanism was obtained in all 
three crystal structures and is considered in detail for the 
Pbnm case. 

The most interesting phase is Pbnm since both the 
HS-LS and metal-insulator transitions occur simultane¬ 
ously with the Pbnm to Pm3m phase transition. The 
GGA+DMFT spectral functions for different cell vol¬ 
umes are shown in Fig. [2j Note, that the Pbnm phase 
doesn’t exist to 2 GPa at 770 K. But BiFeOa in the 
Pbnm structure and experimental with cell volume of 
R3c at ambient pressure is an insulator with a gap of 1 
eV which is close to the experimental data3£~—. When 
the cell volume is 0.85Vo the gap is closed but iron is 
still in the HS state (see Fig. [3]) that contradicts to the 
mechanism proposed in Ref. [2|. Then together with the 
volume decrease the t 2 g occupancy increases and both 
magnetic moment and a number of e g electrons shrink. 
Both transitions don’t appear to happen instantly but 
have a crossover region 0.85-0.82Vo that agrees with the 
experimental observation s 2 ! 3 ! 14 . 

MIT at room temperature where long range order ex¬ 
ists occurs in the pressure range 45-55 GPa. In order to 
reproduce the temperature dependence of critical pres¬ 
sure value the paramagnetic GGA+DMFT calculation 
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FIG. 3: (Color online) Magnetic moments (red square), oc¬ 
cupancy of the d states (black circles), f 2 3 (blue diamonds), 
and e g (green triangles) states obtained in the GGA+DMFT 
(CT-QMC) calculations for the Pbnm phase. 


for /3=40eV -1 eV was carried out. The results obtained 
show that at 33 GPa BiFeOa is already in the metallic 
LS state as well as at higher temperature. Hence increase 
of the critical pressure isn’t a temperature-driven effect. 
Then we have carried out magnetic GGA+DMFT cal¬ 
culation for the same cell volumes of Pbnm structures. 
Small external magnetic field 0.01 eV was applied to each 
Fe atom so that direction of the field on sites corresponds 
to experimentally observed G-type AFM. The gap value 
increases slightly but more important magnetic solution 
was stabilized and MIT occurs between 0.82-0.8 Vo that 
corresponds to 33-43 GPa. This value is close to tran¬ 
sition pressure 45-55 GPa measured at 300 K in AFM 
phase. 

Probabilities of charge and orbital configurations mea¬ 
sured in the GGA+DMFT method within impurity 
solver are shown in FigJS] The hybridization between 
iron and oxygen and hence the d states occupation in¬ 
creases with the volume contraction. At ambient pres¬ 
sure, the d 5 configuration has the larger impact but the 
d 6 probability is also sizable. With the volume reduction 
probabilities of d 6 and d 7 configurations grow up and the 
probability of the d 5 configuration shrinks, see the upper 
panel of Figj5] 

More detailed picture of orbital configuration impacts 
is shown in Fig. [5] (lower panel). One can see that the 
probabilities of the HS configuration with any number 
of electrons (d 5 , d 6 and d 7 ) which dominate at ambi¬ 
ent pressure goes down and after the cell volume become 
smaller than 0.85Vo drops and LS and IS configuration 
probability increase. 

The calculation carried out for the high temperature 
cubic Pm3m phase demonstrates that for all experimen¬ 
tal volumes BiFeC >3 has a metallic behavior. The high 
spin to low spin state transition was also obtained in the 
calculation for the Pm3m phase at 0.8 Vq (33 GPa). 
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FIG. 4: (Color online) Combined phase diagram of BiFeOs. 
Experimental data taken from from RefQi are shown with 
solid lines. Insulator and metallic solutions obtained in 
GGA+DMFT calculation are shown with black squares and 
blue triangles, correspondingly. Critical pressure obtained 
in paramagnetic calculation is 25 GPa. Spin polarized 
GGA+DMFT results at 300 K results in IMT in a range 33-43 
GPa. 


IV. IN CONCLUSION 

We have performed the GGA+DMFT calculations for 
BiFeC >3 at 700 K where it has no long range order for 
a number of different cell volumes corresponding to the 
whole range of experimentally applied pressures (from 
ambient pressure up to 60 GPa). In agreement with ex¬ 
periments, the spectral functions for BiFeOa in the i?3c 
phase as well as in Pbnm at ambient pressure demon¬ 
strates an energy gap of about 1.2 eV in a good agree¬ 
ment with experiment. When the unit cell volume is 
0.85-0.82 Vo that corresponds to 25-33 GPa BiFe 03 in 
the Pbnm crystal structure becomes metallic, and, si¬ 
multaneously, its magnetic moments shrink. In the low 
pressure R3c structure, the same transition occurs in cell 
with the smaller volume 0.8-0.76 Vo (40-50 GPa). The 
high pressure cubic Pm3m phase is metallic in the whole 
range of the investigated pressures (cell volumes) but the 
HS-LS transition also takes place when the cell volume is 
smaller than 0.82 Vo (33 GPa). The higher critical pres¬ 
sure of 43 GPa was obtained in the AFM GGA+DMFT 
calculations for the Pbnm structure in good agreement 
with the experimental phase diagram. 
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FIG. 5: (Color online)Probabilities of different charge states 
d n (upper panel) and orbital configurations of Fe ion (lower 
panel) vs. pressure obtained in the GGA+DMFT (CT-QMC) 
calculations for Pbnm phase (shown only largest ones): d 5 - 
HS (up: 44- down: t% g e°), d 5 -LS (up: 44- down: t 2 g e° g ), 
d 6 -HSi (up: 44 down: 44)’ d 6 -HS 2 (up: 44’ down: 
44 )’ d 6 -IS ( U P : 44 , down: t\ g e° g ), d 6 -LS (up: ti g e° g , down: 

44 )’ d 7 -HS (uP : 44 down: 44 )’ rf7 -ISi (up: 44 - 
down: 44 ), d 7 -IS 2 (up: 44 , down: 44 )- 


V. ACKNOWLEDGMENTS 

The authors thank P. Werner for the CT-QMC im¬ 
purity solver. The GGA+DMFT calculations were per- grant No. RUP1-7077-EK-12 (12-CD-2). 


1 Y. Tokura, Science 312 , 1481 (2006). 


2 A.G. Gavriliuk et al, Phys. Rev. B 77 , 155112 (2008). 
























5 


3 I.S. Lyubutin, A.G. Gavriluk, V.V. Struzhkin, JETP Lett 
88 524 (2008). 

4 A.G. Gavriluk et al., JETP Lett 86 197 (2007). 

5 M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

6 V.I. Anisimov et al., J. Pliys.: Condens. Matter 9, 7359 
(1997); A. I. Lichtenstein and M. I. Katsnelson, Phys. Rev. 
B 57, 6884 (1998); K. Held et al, Phys. Stat. Sol. (b) 243, 
2599 (2006). 

1 J. Kunes et al., Nature Materials 7, 198 (2008). 

8 J. Kunes et al, Phys. Rev. Lett. 102, 146402 (2009). 

9 A.O. Shorikov et al, Phys. Rev. B 82, 195101 (2010). 

10 V. V. Mazurenko et al, Phys. Rev. B 81, 125131 (2010). 

11 M. Moreau et al, J. Phys. Chem. Solids 32, 1315 (1971). 

12 D. C. Arnold et al, Phys. Rev. Lett. 102, 027602 (2009). 

13 S. A. T. Redfern et al, arXiv:0901.3748, 2009. 

14 G. Catalan, J. Scott, Advanced Materials. 21, 2463-2485 
(2009). 

15 A. G. Gavriliuk et al, JETP Lett.82, 224 (2005). 

16 O. E. Gonzalez-Vazquez, J. Iniguez, Phys Rev. B 79, 
064102 (2009). 

17 S. Baroni, S. de Gironcoli, A. D. Corso, and P. Giannozzi, 
http://www.pwscf.org 

18 G. H. Wannier, Phys. Rev. 52, 191 (1937). 


19 N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 
(1997); W. Ku, H. Rosner, W. E. Pickett, and R. T. Scalet- 
tar, Phys. Rev. Lett. 89, 167204 (2002). 

20 Dm. Korotin et al., Euro. Phys. J. B 65, 1434 (2008). 

21 A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. 
Rev. B 52, R5467 (1995). 

22 P. H. Dederichset al, Phys. Rev. Lett. 53, 2512 (1984); 
O. Gunnarssonet al. Phys. Rev. B 39, 1708 (1989); V. I. 
Anisimov and O. Gunnarsson, ibid. 43, 7570 (1991). 

23 S.A. Gramsch, R.E. Cohen and S.Yu. Savrasov, Amer. 
Miner. 88, 257-261 (2003). 

24 P. Werner et al, Phys. Rev. Lett. 97, 076405 (2006). 

25 Mark Jarrell and J. E. Gubernatis, Phys. Rep. 269, 133 
(1996). 

26 P. Ravindran et al Phys. Rev. B 74, 224412 (2006). 

27 B. Ramachandrana, and M. S. Ramachandra Rao, J. Appl. 
Phys. 112, 073516 (2012). 

28 R. V. Pisarev, A. S. Moskvin, A. M. Kalashnikova, and 
Th. Rasing, Phys. Rev. B 79, 235128 (2009). 

29 V. Fruth et al J. Eur. Ceram. Soc. 27, 937 (2007). 

30 R. Palai et al, Phys. Rev. B 77, 014110 (2008). 

31 T. P. Gujar, V. R. Shinde, C. D. Lokhande, Mater. Chem. 
Phys. 103, 142 (2007). 


